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By means of ab initio calculations based on the density functional theory we investigated the 
magnetic phase diagram of ordered FePds alloy as a function of external pressure. Considering 
several magnetic configurations we concluded that the system under pressure has a tendency toward 
noncollinear spin alignment. Analysis of the Heisenberg exchange parameters Jij revealed strong 
dependence of iron-iron magnetic couplings on polarization of Pd atoms. To take into account that 
effect we built an extended Heisenberg model with higher order (biquadratic) terms. Minimizing the 
energy of this Hamiltonian, fully parametrized using the results of ab initio calculations, we found 
a candidate for a ground state of compressed FePd3 , which can be seen as two interpenetrating 
triple-Q phases. 



I. INTRODUCTION 

Magnetism of Fe alloys is a long standing problem in 
solid state physics having fundamental importance for 
modern technology.^ One of the big challenges is under- 
standing magnetism of iron on a face-centered cubic (fee) 
lattice, which is a key issue toward comprehension of a 
variety of phenomena, such as spin glass behavior, the In- 
var anomaly, and a quantitative description of structural 
transition in steels. 

In pure 7-Fe first-principles calculations predict the 
stabilization of a spin spiral state. ^ More generally it 
has been shown that, for a number of disordered fee Fe- 
based binary alloys, such as Fe-Ni and Fe-Pt, the sign of 
the effective magnetic interaction (ferro- or antiferromag- 
netic) is dependent on the volume, being ferromagnetic 
at larger volumes and antiferromagnetic at lower ones.^^ 
The antiferromagnetic (AFM) interactions on frustrated 
and chemically disordered lattices lead to the appearance 
of non-collinear ground states (GSs) at certain region of 
volumes, as predicted in work by van Schilfgaarde et al^ 
Indeed it was shown theoretically that this volume region 
can be reached under applied pressure. 

Thanks to the use of diamond anvil cells at synchrotron 
facilities it has become possible to measure x-ray spec- 
tra and investigate properties of matter under extremely 
high pressure. In today's experiments the values of ap- 
plied pressure exceed 100 GPa range, giving an oppor- 
tunity to study ordinary compounds in unconventional 
conditions. It is of particular interest to investigate Fe- 
based transition metals alloys under pressure due to the 
above mentioned volume-dependent magnetic peculiari- 
ties of their behaviour. A couple of experiments at ultra- 
high pressures have been performed on disordered Invar 
Fe-Ni alloy s^^^ as well as on FeaC cementite^^. In partic- 
ulary, the results of these experiments suggest the stabi- 
lization of a spin-glass state under pressure in Fe-Ni and 



Fe-Pt alloys^^ and an abrupt change of the magnetic 
state under some pressure ( magnet ovolume instability), 
which is in general consistent with the prediction of a 
sharp variation of exchange couplings with pressure in 
these systems j^i^ii^ 

FePda alloy has recently attracted a strong interest 
due to its Invar behavior. It was observed experimentally 
that under an applied pressure of 7 GPa the system shows 
an anomalously small thermal expansion. Hence the 
considered compound demonstrates a pressure-induced 
Invar effect, reminiscent that of FesoNiso disordered 
alloy.^ Moreover, as it will be shown below, non-collinear 
spin states are stabilized in FePds at low volumes, simi- 
larly to the FesoNiso case. 

Winterrose et al}^ investigated FePds under pressure 
from both experimental and theoretical viewpoints. Re- 
sults of the X-ray diffraction measurements implied that 
under the pressure of 12 GPa the system undergoes a 
significant volume collapse while preserving its CusAu 
(LI2) structure. Moreover, under the same applied pres- 
sure they observed a disappearance of quantum beats in 
a nuclear forward scattering (NFS) experiment, which 
indicates the loss of long-range magnetic order in the 
system. In order to interpret the results, authors per- 
formed a set of supercell calculations, based on Density 
Functional Theory (DFT), for few possible magnetic con- 
figurations. By comparing elastic properties of different 
states with experimental data, they came to the con- 
clusion that there is a high-spin (HS) to low-spin (LS) 
transition taking place under pressure. This hypothesis 
is supported by the fact that the obtained magnetic mo- 
ments in the LS state are of the order of 10~^ jiB and 
this can explain the loss of signal in NFS experiment. On 
the other hand, the LS phase never had lower total en- 
ergy than the ferromagnetic (FM) state in a considered 
volume range. However, the possibility of paramagnetic 
case, i.e. the Curie point under pressure being below 
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room temperature, cannot be simply excluded. 

FePda is a complex magnetic system, because it is com- 
posed of rather localized (Fe) and itinerant (Pd) magnetic 
moments. One of the first attempts to account for coex- 
isting magnetism of different characters was carried out 
by Mohn and Schwarz.^^ They proposed a model where 
local spins produced an effective Weiss field acting on 
itinerant magnetic sublattice. The parametrization of the 
model was based on the results of ab initio calculations. 
The developed model was applied for Pd-rich Fe^^Pdi-a^ 
{x < 0.1) alloys and estimated Curie temperatures (Tc) 
were found to be in good agreement with experiment. 

The problem of induced moments and their influence 
on magnetic properties of some Fe-containing alloys was 
studied in a series of publications by Mryasov et ali^^ 
For example, it was pointed out that polarized Pt atoms 
are responsible for anomalous temperature behavior of 
magnetocrystalline anisotropy in FePt. And the induced 
magnetization on Rh atoms was shown to play a crucial 
role in the phase stability of FeRh. 

An extensive first-principles study of iron-palladium 
compounds was made by Burzo et al.^^ There authors 
performed calculations using scalar relativistic tight- 
binding linear muffin-tin orbital method (TB-LMTO)^^ 
The reported values of equilibrium lattice constant and 
magnetic moments were in good agreement with avail- 
able experimental data. A rough estimation of Tc for 
different Pd concentrations using mean field approxima- 
tion (MFA) was not successful, yielding a too low value 
(while MFA should normally overestimate the result). 

A more accurate approach for calculating Tc in Fe-Pd 
solutions and compounds was proposed in Ref. [2l|. Pd 
magnetization was considered to be proportional to the 
vector sum of neighboring magnetic moments. The sug- 
gested computational scheme was based on the extended 
Heisenberg Hamiltonian with exchange parameters ex- 
tracted from self-consistent ab initio calculation. 

Recently another group reported a theoretical study of 
the magnetism of FePda under applied pressure^^ The 
authors explored various chemically and magnetically 
disordered states and came to the conclusion that these 
states can not be the candidates for a GS of the system. 
On the other hand they observed a strong competition 
between commensurate (FM,AFM) and incommensurate 
[spin spiral (SS)] magnetic configurations. The authors 
suggest that this is an indication that the system might 
undergo a transition to some noncollinear state upon the 
application of pressure. 

Mainly motivated by the results of Winterrose et a/.,^^ 
in the present work we will show a possible GS of or- 
dered FePda under pressure and discuss the origins of 
the magnetic transition. 



II. METHODS 

In order to have a realistic description of electronic 
structure of studied system we carried out a set of self- 




FIG. 1. Schematic representation of ordered FePds. The alloy 
has fcc-based LI2 structure, where Fe atoms are located at the 
corners of the cube and Pd atoms are at the centers of the 
faces. 



consistent DFT calculations using PY-LMTO code^^ 
Correlation effects were treated within Local Spin Den- 
sity Approximation (LSDA) with the parametrization of 
Vosko, Wilk and Nusair.^^ 

The crystalline structure of the compound under con- 
sideration is depicted in Fig. 1. In the present study we 
were interested in finding a magnetic GS of a compressed 
FePda. It is certainly impossible to explore all magnetic 
configurations as there are an infinite number of them. 
Thus we have investigated the most plausible candidates 
for the GS by comparing their total energies and examin- 
ing their stability with respect to infinitesimal deviations 
of magnetic moments orientations, keeping in mind that 
this is a probe of the local stability and not the global 
one. 

First we considered possible SS states. These simu- 
lations were carried out on the basis of the generalized 
Bloch theorem^ allowing to avoid laborious supercell 
calculations. Note, that these calculations are scalar 
relativistic and therefore there is no coupling between 
spin degrees of freedom and the crystal lattice. Due 
to the presence of a global spin rotational invar iance, 
the direction of quantization axis can be chosen arbitrar- 
ily. We will align it along the z axis. Hence SS states 
are defined by four parameters: the propagation vector 
[Q = {Qx^Qv^Qz)] and the angle 9 formed by magnetiza- 
tion and the z axis. Once these parameters are chosen, 
the magnetization of the iron atom in neighbouring cells 
is rotated by an angle = Q • R, where R denotes a 
translational vector. Initial phases of Pd moments lo- 
cated at positions were set to (/>^=q • t^, but were 
allowed to choose preferred orientation during the self- 
consistent loop. The SS calculations are shown in Sec. 
IV. 

Consideration of more complex magnetic phases re- 
quired construction of appropriate super cells. In order to 
accommodate some spin structures, we had to use cells 
containing up to eight formula units. The directions of 
iron magnetic moments were prescribed and frozen dur- 
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ing these calculations. This was done so because the 
differences in energies associated to spin deviations are 
rather small (~meV), so the total energy profile is very 
shallow in this direction and in addition possesses plenty 
of local minima. Due to induced character of Pd mag- 
netic moments, their directions and magnitudes were ob- 
tained fully self-consistently. The corresponding calcula- 
tions are described in sections V and VI. 

The high-temperature paramagnetic (PM) state was 
modeled by a disordered local moment (DLM)^ configu- 
ration. The effect of magnetic disorder was taken into ac- 
count by the coherent potential approximation (CPA), as 
implemented in TB-LMTO-CPA.^^ The DLM approach 
is used to describe properties of the system consisting of 
randomly distributed magnetic impurities embedded in a 
non-magnetic medium. Thus induced magnetization on 
a non- magnetic ions, which is a result of an overlap with 
spin-polarized bands originating from magnetic atoms, 
collapses to zero in this phase. Therefore, such an ap- 
proach can be a useful tool for understanding the nature 
of magnetic moments. 

In order to get a deeper insight into magnetic prop- 
erties, we mapped our system on a classical Heisenberg 
model of the following form: 
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FIG. 2. The Fe-Fe exchange interactions in FePds as functions 
of the relative inter-atomic distance d/a, where a is the lattice 
parameter. The results obtained from the DLM reference 
state calculated for different cell compressions. 



few per cent. To avoid any ambiguity, we will work with 
the experimental lattice parameter for FePda, which is 
3.849 A at ambient conditions, ^^^^^ and its correspond- 
ing volume Vo is used throughout the paper. 



where Jij denotes the exchange integral between mag- 
netic atoms at sites i and j, and and ej are unit vec- 
tors in the directions of the local magnetization on sites 
i and j, respectively. Exchange parameters were com- 
puted using the approach of Lichtenstein et al based on 
magnetic force theorem: 

Jij = /_2 TriA,GlA,G],)de (2) 

where is an exchange potential on the i-th site and G^j 
is an intersite Green's function, which describes propaga- 
tion of an electron with spin a = {t^i} from site i to j. 
Within this method, exchange integrals between two sites 
are calculated as a response on infinitesimally small de- 
viations of corresponding magnetic moments away from 
the reference state. So the assumed magnetic order mat- 
ters and, as it will be shown later, the effective Jij pa- 
rameters can be different for various states. The conse- 
quences of this will be discussed in more details in Sec. 
III. In order to probe the stability of a certain magnetic 
state, extracted J's were used to find low-energy mag- 
netic excitation (i.e. frozen magnon) spectra, which are 
the eigenmodes of the considered Hamiltonian. More de- 
tailed description of computational aspects can be found 
elsewhere. 

It has to be mentioned that equilibrium volume value 
and corresponding pressures predicted by LSDA are dif- 
ferent from the experimental ones. In most of the cases 
LSDA leads to underestimation of the bond length by a 



III. HEISENBERG EXCHANGE 
INTERACTIONS 

Self-consistent calculations for the DLM state were car- 
ried out in TB-LMTO-CPA. It is found that Fe atoms 
keep the values of their moments in DLM(PM) state al- 
most unchanged as compared to the magnetically ordered 
states, thus suggesting a high degree of localization of the 
iron magnetic moment. The difference in absolute values 
of Mi?e in DLM and FM states at the equilibrium volume 
was approximately 2%. The same correspondence holds 
under applied pressure and thus there is no tendency to- 
ward a drop of the magnetic moment as suggested by the 
HS-LS scenario. Similar results were already reported 
in a previous study^ 

Using obtained TB parameters, we have calculated 
pair exchange integrals starting from the DLM state. In 
this case Fe moments do not have any prescribed ori- 
entation and extracted J— parameters should reflect the 
properties of the system in the high-temperature phase 
above Tc. This approach was successfully applied in pre- 
vious studies j^i^ 

It should be mentioned, however, that application of 
the magnetic force theorem to systems with induced lo- 
cal moments has limited validity, as was shown by San- 
dratskii et al.^^ In the DLM state small induced Pd mo- 
ments are reduced to zero due to random orientations of 
Fe spins. We therefore investigate a net effect of iron 
moments on possible magnetic ground states which may 
exist at low temperatures. 

In Fig. [2] we show calculated effective JpeFe parame- 
ters as a function of inter-atomic distance for a few frac- 
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FIG. 3. Iron- iron exchange parameters obtained from DLM 
and AFM[110] states for O.SVq. Two iron sublattices which 
appear in the AFM[110] are denoted as Fca and Fes. 



tions of equilibrium volume. The results suggest that 
the first and second nearest-neighbor (NN) interactions 
are dominant. While the first NN interactions (6 neigh- 
bors) are FM, the second NN (12 neighbors) couplings 
are AFM. Third NN exchange parameters, which are FM, 
are also important. Such oscillatory behaviour is due to 
the Ruderman-Kittel-Kasuya-Yosida (RKKY)^— nature 
of exchange interactions in metals and is similar to that 
reported for bcc Fe.^^ 

We report a strong increase of second NN AFM inter- 
actions with pressure while all other couplings depend on 
volume much more weakly. Moreover, the corresponding 
neighborhood forms the fee lattice, which is frustrated for 
this sign of interaction. Here we found that the frustra- 
tion, being a natural source of noncollinearity in spin sys- 
tems, effectively increases its contribution at lower vol- 
umes. We will refer to this fact in the next section. 

Next we have examined the stability of several mag- 
netic states in a compressed FePda. In Fig. [3] we show 
the calculated iron-iron JpeFe exchange parameters ex- 
tracted from DLM and AFM [110] states. Using obtained 
parameters, we computed frozen magnon spectra for each 
of these states. However, all phases demonstrated lo- 
cal magnetic instabilities, indicating that any small per- 
turbation would destroy the state. In other words all 
these magnetic configurations are not even local minima 
on a phase diagram of the system under consideration. 
Nonetheless, we obtained important information: the Fe 
spins have a tendency toward noncollinearity. Motivated 
by this fact, we examined the energies of spin spiral and 
other noncollinear states. 



IV. SPIN SPIRAL CALCULATION 

First we explored the manifold of the states charac- 
terized by O = 90^. The energies of these states as a 
function of a wave vector are shown in Fig. HJ The re- 
sults imply that for ambient pressure F point has the 
lowest energy among all considered configurations. This 



is a manifestation of the stability of the FM state. 
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FIG. 4. Total energies of spin spiral states relative to the 
energy of the FM state at a given volume in FePds. R—X and 
R — M directions are not shown, as wave- vector dependence of 
the total energy was found to be monotonic along these paths. 
The energies are given per chemical formula unit (FU). 

Nevertheless, there are two local minima at X and 
M high symmetry points. These points correspond to 
antiferromagnetic states with ordering vectors [100] and 
[110], respectively. One additional minimum is lying 
along T — R direction and corresponds to wave vector 
q=(^, ^)- Such peculiarities of the total energy 
were already emphasized by another research group. 
The overall shape of the energy profile is in fair agree- 
ment with those calculations, even though different com- 
putational schemes were used. 

It is seen that with increasing pressure the stability of 
the FM state is reduced and at the volume about O.88V0 
we observe the magnetic transition at the point M cor- 
responding to the AFM [110] phase, as was previously re- 
ported by Winterrose et al}^ Further compression leads 
to further destabilization of the ferromagnetic configura- 
tion. 

It should be noted that iron has a quite rigid magnetic 
moment in the entire considered volume range: For ma- 
jority of configurations, as volume is decreased by 20%, 
its magnetization Mp^ lowers by not more than 11%. 

Another remarkable fact is that, for a fixed volume, the 
magnitude of the iron magnetic moment has very similar 
values among different SS states. The highest difference 
was observed between the FM and AFM [110] states and 
was estimated to be ^0.1 (Ib per atom. Meanwhile the 
value of Mpd strongly depends on its environment as 
was already pointed out by another group. For exam- 
ple, at Vo volume in FM state all Pd ions possess the 
magnetic moment of ^0.35 jiB per atom. In the lay- 
ered AFM [100] state Pd atoms which belong to the same 
layer as Fe atoms have a magnetization of 0.14 jiB per 
atom, pointing parallel to the iron moment. The remain- 
ing palladium atoms are non-magnetic. In the AFM [110] 
all Pd moments collapse to zero, because each of them is 
surrounded by an equal number of Fe moments pointing 
"up" {Fca) and "down" (Fe^). 
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FIG. 5. Energies of helical spin configurations as a function 
of the angle between magnetization and the z axis in FePds. 
Values are relative to the energy of the FM state at a given 
volume. 
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FIG. 6. Results of the total energy DFT calculations of canted 
states in ordered FePds- The inset shows the magnetic struc- 
ture under consideration. Red arrows indicate iron magnetic 
moments. Pd magnetization is not shown. 



Our next step was an investigation of various helical 
structures. In this set of calculations we have chosen 
two Q-vectors, corresponding to the lowest states, ob- 
served so far, namely AFM[100] and AFM[110]. Freez- 
ing these three parameters of the SS, we tried to vary the 
value of the O angle. Energies of such magnetic struc- 
tures are shown in Fig. [5l The results again suggest a 
reduction of the FM stability with increasing pressure, 
but for both studied directions (Q), we observed a wide 
range of volumes where helical states are in favor. As 
one compresses the cell up to 0.92Vo the FM phase be- 
comes almost degenerate with two more states, which in 
the coordinate system (Q, 9) correspond to the points 
(X,50^) and (M,50^). Further volume decrease leads to 
the destabilization of the FM solution and the states be- 
longing to the family (M,a) possess the lowest energy in 
the entire volume range. 

At 0.8Vb the AFM[110] phase possesses the lowest en- 
ergy among spin-polarized states. However, we have al- 
ready pointed out that the calculated frozen magnon dis- 
persion indicated a local instability of this state, i.e. ap- 
pearance of imaginary eigenvalues in the corresponding 
spectrum. Therefore we deduce that this is not a OS and 
this fact forced us to investigate more complex states. 
Moreover, so far we have considered only coplanar spin 
structures. 



V. CANTED SPIN STATES 

In order to explore more complex magnetic states, we 
have constructed a 2x2x1 supercell of FePda with four 
uneqivalent iron sublattices. Here we introduce an angle 
d for iron spins, defined in a way shown in the inset of 
Figure [6l Thus ^ = corresponds to the AFM[110] phase 
and Q = 180^ - to AFM[100]. By tuning 6 one can go 



continously from one state to another. At each volume 
the calculations were carried out for the angle 6 fixed to 
a given value. A fully unconstrained determination of 
would be an extremely difficult task, due to the reasons 
mentioned in Sec. II. 

According to the obtained results, corresponding func- 
tion Etot(^) in addition to the points = {0, 180^} has 
one more minimum at certain angle (^o)- Such angu- 
lar dependence can not be accounted within classical 
Heisenberg picture, which should give a "cos (^)" cur- 
vature. Note, that the present shape of the total en- 
ergy profile was found even without inclusion of spin- 
orbit coupling (SOC). Thus relativistic effects, such as 
magnetic anisotropy, are not responsible for such behav- 
ior. The situation is reminiscent of the study of iron 
pnictides,^^ where similar constrained calculations were 
performed. For the latter case it was proposed that the 
proper Etot{0) dependence can be obtained by introduc- 
ing a higher-order exchange terms to the spin model. 

Following this idea, we have built an extended 
Heisenberg-like model, which aims to catch the essential 
physical properties of the FePds and guide us toward 
finding a true ground state. 



VI. EXTENDED HEISENBERG MODEL AND 
TRIPLE-Q STRUCTURE 

The influence of Pd moments on stability of the FM 
state is crucial as was demonstrated by Polesya et alS^ 
and therefore it is essential to take them into account. 

Palladium is a rather peculiar element: it is known 
that already bulk Pd is characterized by a high density 
of states at the Fermi level, so the Stoner criterion is 
nearly satisfied. Therefore Pd is easily polarized by the 
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contact with neighboring magnetic moments. Certainly 
it participates in magnetic interactions and as a result 
effective couplings between two iron sites in FM phase, 
where Pd gets polarized, and in AFM one, where it is 
nonmagnetic, are considerably different (Fig. [3j). 

The underlying physics can be already understood by 
looking on the expression of the J-parameters in Eq. ([2]). 
Since the magnitude of Fe magnetic moment is almost 
configuration-independent, one gets that lA^^el is approx- 
imately the same in all states. Therefore what gives rise 
to the observed difference in exchange couplings is the 
intersite Green's function. In the FM state an electron, 
going from one iron site to another, propagates through a 
strongly polarized medium, while in the AFM state this 
polarization is missing. Since Pd magnetization is large 
(~0.3 jj^B per atom), being a first-order term in A^e, 
its disappearance gives a significant impact on Gij and 
eventually on exchange integrals. As a result, a classical 
Heisenberg model with pairwise interactions and bilinear 
exchange only, cannot properly map the dependence of 
the total energy on the magnetic configurations and must 
be extended to include higher-order exchange terms. 

Such behaviour was already reported for FeRh 
compoundJ^ It is worth emphasizing that this situation 
is not generic and can be ascribed to be the feature of 
these 4(i elements. In order to stress this point, we per- 
formed additional calculations for a hypothetical system, 
where palladium was substituted by copper atoms within 
same geometry of the unit cell. Opposite to the pre- 
vious case, Cu is non-magnetic in both AFM [110] and 
FM states. In this case we found that J2 parameters 
extracted from both configurations almost coincide with 
each other. Thus we do not observe any pronounced de- 
viations from the Heisenberg magnet behaviour. 

Hence, in order to take into account such peculiari- 
ties of the magnetic interactions, we propose an effective 
model for iron degrees of freedom, which in addition to 
Eq. ([1]) contains higher-order exchange terms, originating 
from the polarization of palladium atoms: 

H = Hexch-^Jij{ei-ejf (3) 

where J' is a biquadratic exchange parameter. 

Thus we are able to build a Heisenberg-like model and 
fully parametrize it using the results of DFT calculations. 
For simplicity we considered interactions with first three 
coordination spheres only (Ji, J2,-^3) as the remainder 
are much smaller. As was already pointed out, the Pd- 
originated renormalization of J-parameters is the most 
pronounced for next NNs (J2) and therefore we intro- 
duce biquadratic term only for this coupling. The val- 
ues which were extracted from the Figl3] and used for 
the model are the following: { Ji, J2, ^3} ~^ {0.236,- 
0.902,-0.282,0.29} (in mRy). It is seen that is of the 
same order as bilinear exchange. Such situation is atyp- 
ical, but not unique: for example, a sizable value of the 
biquadratic term is necessary to explain properties of an- 
other class of Fe-based mater ials.^^ Moreover, the present 




FIG. 7. Obtained high-pressure magnetic phase of FePds. 
Arrows represent iron magnetic moments. Four spins of the 
same color, forming the triple- Q state, point toward the center 
of the cube, four others - opposite to it. Note, that the state is 
degenerate with respect to simultaneous rotations of all spins 
through the same angles. 



set of parameters describes well the curvature of Etot{0) 
profile, shown in Fig. 6. The very existence of the mini- 
mum ^0 as well as its position are in fair agreement with 
those calculations. This can be viewed as an indication 
of plausability of the chosen parameters for our model. 

Minimizing the energy of the Hamiltonian (Eq. (3)) on 
a 2x2x2 supercell we obtained a new ground state of 
the system, which is depicted in Fig. [3 The magnetic 
structure of this state can be seen as two interpenetrating 
fee subsystems, whose spins form so-called triple-Q (3Q) 
states. The angle between each pair of spins is 109^28^ 
hence the vectors point toward the vertices of an ideal 
tetrahedron. Each spin out of eight has an identical pair, 
so its vertices are doubly degenerate. Analysis of the 
spin arrangement revealed that J3 coupling (FM), which 
connects two fee sublattices, is fully satisfied within such 
geometry. The rest of the interactions compete with each 
other due to presence of frustration in the system. Note 
that negative sign of J2 means that the interaction favors 
perpendicular spin orientation. This term is the driving 
force lifting the degeneracy associated with the frustra- 
tion, and results into stabilization of a certain angle be- 
tween spins. 

It was shown before for fcc-based alloys that stability of 
triple-Q state with respect to single-Q ones can be related 
to nonlinear spin interactions in the system^^ and/or the 
presence of paramagentic impurities^^ . In fact, both of 
these ingredients are present in FePda, and hence our 
findings are consistent with the established physical pic- 
ture. However, to the best of our knowledge the pre- 
sented spin structure was not observed in alloys with such 
a low concentration of magnetic ions. 

In order to check the stability of the obtained config- 
uration, we have computed the corresponding dispersion 
of low-energy magnetic excitations (FiglH]). First of all 
we confirm that all excitations have a positive energy, 
meaning that the state is stable. It is worth emphasiz- 
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FIG. 8. Frozen magnon energies of the proposed high-pressure 
magnetic phase of FePds- 



ing once again, that all states which were studied before 
showed instabilities in their spectra. Second, one can see 
few almost degenerate Goldstone modes obeying linear 
dispersion law at small values of the wave vector q. We 
will refer to this fact in the next section. 

Finally, we carried out total energy DFT calculation of 
the 3Q state for different fractions of equilibrium volume. 
It was confirmed that this noncollinear state becomes 
lower in energy than FM one at a relative compression 
rate of 0.96, which is in excellent agreement with the 
experiment. The triple-Q state possesses the lowest en- 
ergy among all studied states in the low- volume region4^ 

Hence we see that the model with the present choice of 
J— parameters describes well the properties of the system 
next to the OS. It has to be mentioned, however, that this 
set is unable to reproduce the relative energies of spin 
spirals depicted in Fig. HI One of possible explanations 
is that incommensurate spin states possess a significant 
Pd magnetization which has to be taken into account 
explicitly for a proper energy estimation. 



VII. OBSERVATION OF THE TRIPLE-Q STATE 

3Q states are rather difficult to observe experimen- 
tally, because their behavior is similar to collinear 
antiferromagnets.^^ The measurement which helped to 
distinguish these two phases was proposed by Kawarazaki 
et al.^^ The method, however, requires certain elements 
which are the sources of the 7-rays. Thus for the case of 
FePda a more useful way would be to use the Mossbauer 
effect in Fe, but this technique can not provide an un- 
equivocal answer if the state is more complex than single- 

As was already pointed out, first of all we suspect that 
there should be an abrupt change in the shape of the spin- 
wave dispersion from parabolic to linear in FePda under 
applied pressure. This would be the first indication of 



existence of the triple-Q state, which is necessary, but 
not sufficient. 

Another external parameter which has to be controlled 
during an experiment is the temperature. In Ref. |15| 
a complete vanishing of Fe average magnetization was 
found in the NFS experiment upon compression. One of 
explanations is that the value of Tc is strongly affected by 
applied pressure and at a given volume it is pushed below 
room temperature. However, a reliable ah initio evalua- 
tion of Tc of FM FePda would be already a difficult task 
and seems not to be solved in previous studies. ^^^^^ Re- 
cent results on estimation of the Tc in MFA confirm its 
decrease upon the application of pressure^^, but it seems 
to reach the room temperature at lower volumes than 
observed experimentally. It should not be excluded that 
3Q configuration might have a different ordering tem- 
perature as compared with the FM one. This is one of 
explanations of the loss of the NFS signal. 

Another issue with such type of systems is disorder. 
Ordered FePda samples are obtained by annealing and 
subsequent fast cooling (quenching). Hence some amount 
of residual disorder can always appear in a real system. 
It was reported that the value of transition pressure dif- 
fers from one sample to another. This fact may suggest 
that a small antisite disorder is present and can affect 
a phase diagram of the compound under consideration. 
In principle, this can give rise to more complicated spin 
structures. 



VIII. CONCLUSIONS 

We have studied theoretically the ground-state mag- 
netic properties of FePda ordered alloy under the external 
pressure using first-principles methods. It was deduced 
that the compound undergoes a magnetic transition from 
the FM to the Triple-Q state at 0.96 fraction of equilib- 
rium volume. Fe atoms possess a significant magnetic 
moment in the high-pressure phase, unlike the HS-LS sce- 
nario suggested earlier. The disappearance of the quan- 
tum beats in the NFS experiments^ may indicate a drop 
of Tc across the transition. Essential ingredients which 
stabilize the 3Q states can be determined: (i) an exis- 
tence of strong magnetically frustrated couplings in the 
system, which is J2 in the present case; and (ii) admix- 
ture of higher-order interactions, favoring perpendicular 
spin alignment. 
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